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^ (57) Abstract: The invention concerns the processing of audio data. The invention is characterized in that it consists in: (a) encoding 

^ signals representing a sound propagated in three-dimensional space and derived from a source located at a first distance (P) from 

^5 a reference point, to obtain a representation of the sound through components expressed in a spherical harmonic base, of origin 

f<| corresponding to said reference point, (b) and applying to said components compensation of a near-field effect through filtering 

Q based on a second distance (R) defining, for sound reproduction, a distance between a reproduction point (HPj) and a point (P) of 
auditory perception where a listener is usually located. 
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(57) Abrege : La pnSsente invention concerne le traitement de donn6es sonores. Selon l'invention : a) on code des signaux repr6- 
sentatifs d'un son se propageant dans l'espace tridimensionnel et issu d'une source situee a une premiere distance (p) d'un point de 
reference, pour obtenir une representation du son par des composantes exprim6es dans une base d'harmoniques sph^riques, d'origine 
correspondant audit point de reference, b) et on applique a ces composantes une compensation d'un effet de champ proche par un 
filtrage qui est fonction d'une seconde distance (R) dSfinissant, pour une restitution du son, une distance entre un point de restitution 
(HPj) et un point (P) de perception auditive ou se place g^neralement un auditeur. 
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Procede de traitement de donnees sonores et dispositif 
d 1 acquisition sonore mettant en oeuvre ce procede 

La presente invention concerne le traitement de donnees 
5 sonores . 

Des techniques relatives a la propagation d*une onde 
sonore dans I'espace tridimensional, impliquant notamment 
une simulation et/ou une restitution sonores specialisees, 

10 met tent en oeuvre des precedes de traitement du signal 
audio appliques a la simulation de ph^nomSnes acoustiques 
et psycho- acoustiques. De tels proc^des de traitement 
prevoient un encodage spatial du champ acoustique, sa 
transmission et sa reproduction spatialisee sur un 

15 ensemble de hauts-parleurs ou sur des ecouteurs d'un 
casque stereophonique . 

Parmi les techniques de son spatialis^, on distingue deux 
categories de traitements compl€mentaires l'une de 1' autre 
20 mais qui sont general ement mise en oeuvre, l'une et 
l 1 autre, au sein d'un meme. systeme . 

D'une part, une premiere categorie de traitements concerne 
les procedes de synthese d'effet de salle, ou plus 

25 generalement d'effets environnementaux . A partir d'une 
description d , une ou plusieurs sources sonores (signal 
emis, position, orientation, directivity, ou autre) et en 
se basant sur un modele d ! effet de salle (impliquant une 
geometrie de salle, ou encore une perception acoustique 

30 souhaitee) , on calcule et l'on dScrit un ensemble de 
phenomenes acoustiques el^mentaires (ondes directes, 
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refl^chies ou dif fractees) , ou encore un phenomene 
acoustique macroscopique (champ r£verb§r€ et dif f us) , 
permettant de traduire l'effet spatial au niveau d'un 
auditeur situ4 a un point choisi de perception auditive, 
5 dans I'espace tridimensional . On calcule alors un 
ensemble de signaux associes typi quern en t aux reflexions 
(sources " secondaires" , actives par re-emission d'une onde 
principale re<?ue, ayant un attribut de position spatiale) 
et/ou associes £ une reverberation tardive (signaux 
10 decorrel^s pour un champ dif f us) . 

D f autre part, une seconde categorie de procedes concerne 
le rendu positionnel ou directiormel • de sources sonores. 
Ces procedes sont appliques a des signaux determines par 
15 un procede de la premiere categorie decrite ci-avant 
(impliquant des sources primaires et secondaires) en 
fonction de la description spatiale (position de la 
source) qui leur est associee. En particulier, de tels 
procedes selon cette seconde categorie permettent 
20 d'obtenir des signaux a diffuser sur des hauts-parleurs ou 
ecouteurs, pour finalement donner a un auditeur 
1' impression auditive de sources sonores plac^es a des 
positions respectives predetermines , autour de 
1' auditeur. Les procedes selon cette seconde categorie 
25 sont qualifies de "createurs d' images sonores 
tridimensionnelles", du fait de la repartition dans 
I'espace tridimensionnel du ressenti de la position des 
sources par un auditeur. Des procedes selon la seconde 
categorie comportent generalement une premiere etape 
d'encodage spatial des evenements acoustiques eiementaires 
qui produit une representation du champ sonore dans 
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l'espace tridimensionnel. Dans une seconde etape, cette 
representation est transmise ou stockee pour un usage 
differe. Dans une troisieme etape, de d£codage, les 
signaux decodes sont deiivres sur des hauts-parleurs ou 
des Scouteurs d'un dispositif de restitution. 

La presente invention s 1 inscrit plutot dans la seconde 
cat£gorie precitee. Elle conceme en particulier 
1' encodage spatial de sources sonores et une specification 
de la representation sonore tridimensionnelle de ces 
sources, Elle s' applique aussi bien a un encodage de 
sources sonores " virtuelles" (applications o\i des sources 
sonores sont simulees telles que des jeux, une conference 
spatialisee, ou autres) f qu ! un encodage "acoustigue" d'un 
champ sonore naturel, lors d'une prise de son par un ou 
plusieurs reseaux tridimensionnel s de microphones. 

Parmi les techniques envisageables de spatialisation du 
son, 1 1 approche "a/nbisonigue " est pr£f £r<§e . L ' encodage 
ambisonique, qui sera d6crit en detail plus loin, consiste 
a reprSsenter des signaux relatifs a une ou plusieurs 
ondes sonores dans une base d'harmoniques spheriques (en 
coordonnees spheriques impliquant notamment un angle 
d' elevation et un angle azimutal, caracterisant une 
direction du ou des sons) . Les composantes representant 
ces signaux et exprimees dans cette base d'harmoniques 
spheriques sont aussi fonction, pour les ondes emises en 
champ proche, d'une distance entre la source sonore 
•emettant ce champ et un point correspondant & l'origine de 
la base des harmonique spheriques. Plus particulierement , 
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cette dependance de la distance s'exprime en fonction de 
la frequence sonore, comme on le verra plus loin. 

Cette approche ambisonique offre un grand nombre de 
5 fonctionnalites possibles., notamment en terme de 
simulation de sources virtuelles, et, de manidre generale, 
pr^sente les avantages suivants : 

- elle traduit, de fagon rationnelle, la rgalite des 
phenomenes acoustiques et apporte un rendu audit if spatial 

10 realiste, convaincant et immersif ; 

- la representation des phenomdnes acoustiques est 
scalable : elle offre une resolution spatiale qui peut 
§tre adaptee k diff6rentes situations. En effet, cette 
representation peut §tre transmise et exploits en 

15 fonction de contraintes de d£bit lors de la transmission 
des signaux encodes et/ou de limitations du dispositif de 
restitution ; 

- la representation ambisonique est flexible et il est 
possible simuler une rotation du champ sonore, ou encore, 

20 a la restitution, d' adapter le decodage des signaux 
ambisoniques a tout dispositif de restitution, de 
geometries di verses. 

Dans 1' approche ambisonique connue, l'encodage des sources 
15 virtuelles est essentiellement directionnel . Les fonctions 
d'encodage reviennent a calculer des gains qui dependent 
de l 1 incidence de l'onde sonore exprimee par les fonctions 
harmoniques sph£riques qui dependent de I 1 angle 
d' elevation et de 1' angle azimutal en coordonnees 
10 spheriques. En particulier, au decodage, on suppose que 
les hauts-parleurs, a la restitution, sont lointains. II 
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en resulte une distorsion (ou line incurvation) de la forme 
des fronts d'onde reconstruits . En effet, comme indiqug 
ci-avant, les composantes du signal sonore dans la base 
des harmoniques spheriques, pour un champ proche, 
5 dependent en fait aussi de la distance de la source et de 
la frequence sonore. Plus precis€ment, ces composantes 
peuvent s 1 exprimer mathematiquement sous la forme d 1 un 
polynome dont la variable est inversement proportionnelle 
k la distance precit^e et 3. la frequence sonore. Ainsi, 

10 les composantes ambisoniques, au sens de leur expression 
theorique, sont divergentes dans les basses frequences et, 
en particulier, tendent vers I'infini quand la frequence 
sonore d^croit vers zero, lorsqu'elles representent un son 
en champ proche emis par une source situSe a une distance 

15 finie. Ce phenomene math^matique est connu, dans le 
domaine de la representation ambisonique, d6ja pour 
l'ordre 1, par le terme de "bass boost", notamment par : 
- M.A.GERZON, "General Met a theory of Auditory 

Localisation", preprint 3306 of the 92 nd AES Convention, 

20 1992, page 52. 

Ce phenomene devient particulierement critique pour des 
ordres d 1 harmoniques spheriques elevees impli quant des 
polynome s de puissance elevee. 

25 On connait par : 

SONTACCHI et HOLDRICH, "Further Jnvestigrations on 3D Sound 
Fields using- Distance Coding 1 ' (Proceedings of the COST G-6 
Conference on Digital Audio Effects (DAFX-01) , Limerick, 
Irlande, 6-8 Decembre 2001) , 
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une technique pour prendre en compte une incurvation des 
fronts d'ondes au sein d'une representation proche d'une 
representation ambisonique, dont le principe consiste a : 
~ appliquer un encodage ambisonique (d'ordre eleve) aux 
signaux issus d'une prise de son virtuelle (simul^e) , de 
type WFS (pour "Wave Field Synthesis") ; 

- et reconstruire le champ acoustique sur une zone d'apres 
ses valeurs sur une frontiere de zone, se fondant ainsi 
sur le principe de HUYGENS - FRESNEL . 



Cependant, la technique presentee dans ce document, bien 
que prometteuse du fait qu'elle utilise une representation 
ambisonique a un ordre eleve, pose un certain nombre de 
problemes : 

- les ressources inf ormatiques n€cessaires pour le calcul 
de toutes les surfaces permettant d' appliquer le principe 
de HUYGENS -FRESNEL, ainsi que les temps de calcul 
necessaires, sont excessifs ; 

- des artefacts de traitement dits d 1 "aliasing spatial" 
apparaissent a cause de la distance entre les microphones, 
a moins de choisir un maillage de microphone virtuels 
serr6 dans l'espace, ce qui alourdit les traitements ; 

. - cette technique est dif f icilement transposable a un cas 
reel de capteurs a disposer en rgseau, en presence d'une 
25 source r<§elle, a 1 1 acquisition ; 

" a la restitution, la representation sonore 
tridimensionnelle est implicitement assujettie a un rayon 
fige du dispositif de restitution car le decodage 
ambisonique doit se faire, ici, sur un reseau de hauts- 
parleurs de m§mes dimensions que le reseau de microphones 
initial, ce document ne proposant aucun moyen d 1 adapter 
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1'encodage ou le decodage aL d'autres tailles de 
dispositifs de restitution. 

Surtout, ce document presente un reseau horizontal de 
capteurs, ce qui suppose que les phenomenes acoustiques 
dont on tient compte, ici, ne se propagent que dans des 
directions horizontales , ce qui exclut toute autre 
direction de propagation et qui, done, ne reprisente pas 
la reality physique d ! un champ acoustique ordinaire, 

De fagon plus g<§nerale, les techniques actuelles ne 
permettent pas de traiter de fagon satisf aisante tout type 
de sources sonores, notamment en champ proche, mais plutdt 
des sources sonores lointaines (ondes planes) , ce qui 
correspond a une situation restrictive et artificielle 
dans nombreuses applications. 

Un objet de la presente invention est de fournir un 
proced€ pour traiter, par encodage, transmission et 
20 restitution, un type quelconque de champ sonore, en 
particulier l'effet d'une. source sonore en champ proche. 

Un autre objet de la presente invention est de fournir tin 
proc6d€ permettant 1' encodage de sources virtuelles, non 
25 seulement en direction, mais aussi en distance, et de 
dSfinir un decodage adaptable a un dispositif de 
restitution quelconque. 

Un autre objet de la presente invention est de fournir un 
30 proc£d£ de traitement robust e pour des sons de toutes 
frequences sonores (y compris les basses frequences) , 
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notamment pour la prise de son de champs acoustiques 
naturels a l'aide de reseaux tridimensionals de 
microphones . 



10 



15 



20 



A cet effet, la pr^sente invention propose un procede de 
traitement de donnees sonores, dans lequel : 

a) on code des signaux representatif s d'au moins un son se 
propageant dans l'espace tridimensional et issu d'une 
source situee a une premiere distance d'un point de 
reference, pour obtenir une representation du son par des 
composantes exprimges dans une base d'harmoniques 
spheriques, d'origine correspondant audit point de 
reference, et 

b) on applique auxdites composantes une compensation d'un 
effet de champ proche par un filtrage qui est fonction 
d'une seconde distance definissSnt sensiblement , pour une 
restitution du son par un dispositif de restitution, une 
distance entre un point de restitution et un point de 
perception auditive. 



Dans un premier mode de realisation, ladite source etant 
lointaine du point de reference, 

- on obtient des composantes d'ordres successifs m pour la 
representation du son dans ladite base d'harmoniques 

25 spheriques, et 

- on applique un filtre dont les coefficients, appliques 
chacun a une composante d'ordre m, s'expriment 
analytiquement sous la forme de 1 1 inverse d'un polynome de 
puissance m, dont la variable est inversement 
proportionnelle k la frequence sonore et a ladite seconde 
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distance, pour compenser un effet de champ, proche au 
niveau du dispositif de restitution. 

Dans un second mode de realisation, ladite source etant 
une source virtuelle prevue a ladite premiere distance, 

- on obtient des composantes d'ordres successifs m pour la 
representation du son dans ladite base d'harmoniques 
spheriques , et 

- on applique un filtre global dont les coefficients, 
appliques chacun a une composante d'ordre m, s'expriment 
analytiquement sous la forme d'une fraction, dont : 

- le numerateur est un polynome de puissance m, dont 
la variable est inversement proportionnelle a la 
frequence sonore- et a ladite premiere distance, pour 
simuler un effet de champ proche de la source 
virtuelle, et 

- le denominateur est un polynome de puissance m, 
dont la variable est inversement proportionnelle a la 
frequence . sonore et a ladite seconde distance, pour 
compenser l 1 effet du champ proche de la source 
virtuelle dans les basses frequences sonores. 

Pref erentiellement , on transmet au dispositif de 
restitution les donnees codees et filtrees aux etapes a) 
et b) avec un parametre representatif de ladite seconde 
distance . 

En complement ou en variante, le dispositif de restitution 
comportant des moyens de lecture d'un support m6moire, on 
memorise sur un support memoire destine a §tre lu par le 
dispositif de restitution les donnees codees et filtrees 
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aux stapes a) et b) avec un parametre representatif de 
ladite seconde distance. 

Avantageusement, prealablement a une restitution sonore 
par un dispositif de restitution comportant une pluralite 
de hauts-parleurs disposes a une troisierne distance dudit 
point de perception auditive, on applique aux donn§es 
codees et filtrees un filtre d' adaptation dont les 
coefficients sont fonction desdites seconde et troisidme 
distances. 

Dans une realisation particuliere, les coefficients de ce 
filtre d' adaptation, appliques chacun a une composante 
d'ordre m, s f expriment analytiquement sous la forme d'une 
fraction, dont : 

- le num^rateur est un polynome de puissance m, dont la 
variable est inversement proportionnelle a la frequence 
sonore et a ladite seconde distance, 

- et le denominateur est un polynome de puissance m, dont 
la variable est inversement proportionnelle a la frequence 
sonore et a ladite troisierne distance. 

Avantageusement, pour la mise en ceuvre de I'gtape b) , on 
prevoit : 

- pour des composantes d'ordre m pair, des filtres 
audionumeriques sous la forme d'une cascade de cellules 
d'ordre deux ; et 

- pour des composantes d'ordre m impair, des filtres 
audionumeriques sous la forme d'une cascade de cellules 

30 d'ordre deux et une cellule suppl€mentaire d'ordre un. 
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Dans cette realisation, les coefficients d'un filtre 
audionumerique , pour une composante d f ordre m, sont 
definis a partir des valeurs numeriques des racines 
5 desdits polyn6mes de puissance m. 

Dans une realisation particuliere, les polynomes precit6s 
sont des polynomes de Bessel. 

10 A 1 1 acquisition des signaux sonores, on prevoit 
avantageusement un microphone comportant un reseau de 
transducteurs acoustiques agences sensiblement sur la 
surface d'une sphere dont le centre correspond 
sensiblement audit point de reference, pour obtenir 

15 lesdits signaux representatif s d'au moins un son se 
propageant dans l'espace tridimensional . 

Dans cette realisation, on applique a l'etape b) un filtre 
global pour, d'une part, compenser un effet de champ 
20 proche en fonction de ladite seconde distance et, d' autre 
part, €galiser les signaux issus des transducteurs pour 
compenser une ponderation de directivite desdits 
transducteurs . 

25 Pr^ferentiellement , on prevoit un nombre de transducteurs 
fonction d'un nombre total choisi de composantes pour 
repr<§senter le son dans ladite base d'harmoniques 
sph6riques . 

30 Selon une caracteristique avantageuse, on choisit & 
l'etape a) un nombre total de composantes dans la base des 
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harmoniques spheriques pour obtenir, a la restitution, une 
region de I'espace autour du point de perception dans 
laquelle la restitution du son est fiddle et dont les 
dimensions sont croissantes avec le norabre total de 
5 composantes . 

Preferentiellement, on pr^voit en outre un dispositif de 
restitution comportant un nombre de haut-parleurs au moins 
egal audit nombre total de composantes. 

10 

En variante, dans le cadre d'une restitution avec synthese 
binaurale ou transaurale : 

- on prSvoit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a vine 

15 distance choisie d'un auditeur, 

- on obtient, pour cet auditeur, une information de 
ressenti attendu de la position dans l'espace de sources 
sonores situees a une distance de reference pr€d€terminee 
de 1' auditeur pour 1 ■ application d'une technique dite de 

20 "synthase binaurale" ou "transaurale" , et 

- on applique la compensation de I'gtape b) avec ladite 
distance de reference sensiblement en tant que seconde 
distance. 

25 Dans une variante o\i 1'on introduit une adaptation au 
dispositif de restitution k deux Scouteurs : 

- on prevoit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a une 
distance choisie d'un auditeur, 

30 - on obtient, pour cet auditeur, une information de 
ressenti de la position dans l'espace de sources sonores 
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situees a une distance de reference predetermine de 
l'auditeur, et 

- prealablement a une restitution sonore par le dispositif 
de restitution, on applique aux donnees codees et filtrees 

5 aux etapes a) et b) un filtre d' adaptation dont les 
coefficients sont fonction de la seconde distance et 
sensiblement de la distance de reference. 

En particulier, dans le cadre d'une restitution avec 
10 synthese binaurale : 

- le dispositif de restitution comporte un casque a deux 
ecouteurs pour les oreilles respectives de l'auditeur, 

- et pr€f erentiellement, s€parement pour chaque ecouteur, 
on applique le codage et le filtrage des Stapes a) et b) 

15 pour des signaux respectifs destines a alimenter chaque 
6couteur, avec, en tant que premiere distance, 
respectivement vine distance separant chaque oreille d'une 
position d'une source a restituer dans l'espace de 
restitution. 

20 

Pref erentiellement, on met en forme, aux Stapes a) et b) , 
un systeme matriciel comportant au moins : 

- une matrice comportant lesdites composantes dans la base 
des harmoniques sph^riques, et 

25 - une matrice diagonale dont les coefficients 
correspondent a des coefficients de filtrage de I'etape 
b), 

et on multiplie lesdites matrices pour obtenir une matrice 
rSsultat de composantes compensees. 

30 



De prSfSrence, a la restitution 
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- le dispositif de restitution comporte une plurality de 
haut-parleurs disposes sensiblement & une mSme distance du 
point de perception auditive, et 

- pour decoder lesdites donnees codees et filtrees aux 
5 etapes a) et b) et former des signaux adaptes pour 

alimenter lesdits haut-parleurs : 

* on forme un systeme matriciel comportant ladite 
matrice r^sultat de composantes compensees et une 
matrice de decodage predetermine , propre au 

10 dispositif de restitution, et 

* on obtient une matrice comportant des coefficients 
representatif s des signaux d 1 alimentation des hauts- 
parleurs par multiplication de la matrice resultat 
par ladite matrice de decodage. 

15 

La presente invention vise aussi un dispositif 
d' acquisition sonore, comportant un microphone muni d'un 
r£seau de transducteurs acoustiques disposes sensiblement 
sur la surface d'une sphere, Selon 1' invention le 
20 dispositif comporte en outre une unite de traitement 
agencee pour : 

- recevoir des signaux emanant chacun d'un transducteur , 

- appliquer auxdits signaux un codage pour obtenir une 
representation du son par des composantes exprimees dans 

25 une base d'harmoniques spheriques, d f origine correspondant 
au centre de ladite sphere, 

- et appliquer auxdites composantes un filtrage qui est 
fonction, d'une part, d'une distance correspondant au 
rayon de la sphere et, d' autre part, d'une distance de 

30 r€f £rence . 
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Preferentiellement, le filtrage effectue par 1' unite de 
traitement consiste, d'une part, a egaliser, en fonction 
du rayon de la sphere, les signaux issus des transducteurs 
pour compenser une ponderation de directivite desdits 
transducteurs et, d 1 autre part, a compenser un effet de 
champ proche en fonction de ladite distance de reference. 

D'autres avantages et caracteristiques de 1' invention 
apparaitront a la lecture de la description d6taill6e ci- 
apres et a l'examen des figures qui 1 ' accompagnent , sur 
lesquelles : 

- la figure l illustre schematiquement un syst^me 
d' acquisition et creation, par simulation de sources 
virtuelles, de signaux sonores, avec encodage, 
transmission, decodage et restitution par un dispositif de 
restitution spatialis€, 

- la figure 2 represente plus precisement un encodage de 
signaux def inis a la fois en intensity et par rapport a la 
position d'une source dont ils sont issus, 

20 - la figure 3 illustre les parametres en jeu dans la 
representation ambisonique, en coordonnees spheriques ; 

- la figure 4 illustre une representation par une metrique 
tridimensionnelle dans un repere de coordonnees 

spheriques, d ' harmoniques spheriques Y%„ de differents 
25 ordres ; 

- la figure 5 est un diagramme des variations du module de 
f onctions radiales j a (kr) , qui sont des f onctions de 
Bessel spheriques, pour des valeurs d'ordre m successives, 
ces fonctions radiales intervenant dans la representation 
ambisonique d'un champ de pression acoustique ; 
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- la figure 6 represente 1 ■ amplification due a l'effet de 
champ proche pour differents ordres successifs m, en 
particulier dans les basses frequences ; 

- la figure 7 represente schematiquement un dispositif de 
restitution comportant une plurality de hauts-parleurs 
HPi , avec le point (reference P) de perception auditive 
precite, la premiere distance precitee (referencee p) et 
la seconde distance precitee (referencee R) ; 

- la figure 8 represente schematiquement les parametres 
mis en jeu dans l'encodage ambisonique, avec un encodage 
directionnel, ainsi qu'un encodage de distance selon 
1 1 invention ; 

- la figure 9 represente des spectres d'energie des 
filtres de compensation et de champ proche simules pour 
une premiere distance d'une source virtuelle p = 1 m et 
une pre -compensation de hauts-parleurs situes a une 
seconde distance R = 1,5 m ; 

- la figure 10 represente des spectres d'energie des 
filtres de compensation et de champ proche simules pour 
une premiere distance de la source virtuelle p = 3 m et 
une pre-compensation de hauts-parleurs situes a une 
distance R = 1,5 m / 

- la figure HA represente une reconstruction du champ 
proche avec compensation, au sens de la presente 

25 invention, pour une onde spherique dans le plan 
horizontal ; 

- la figure 11B, a comparer avec la figure 11A, represente 
le front d'onde initial, issu d'une source S ; 

- la figure 12 represente schematiquement un module de 
filtrage pour adapter les composantes ambisoniques recues 
et pre-compensees a l'encodage pour une distance de 
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reference R en tant que seconde distance, k un dispositif 
de restitution comportant une plurality de hauts-parleurs 
disposes a une troisieme distance R 2 d f un point de 
perception auditive ; 
5 - la figure 13A represente schematiquement la disposition 
d'une source sonore M, a la restitution, pour un auditeur 
utilisant un dispositif de restitution appliquant une 
synthese binaurale, avec une source emettant en champ 
proche ; 

10 - la figure 13B represente schematiquement les Stapes 
d f encodage et de decodage avec effet de champ proche dans 
le cadre de la synthese binaurale de la figure 13A a 
laquelle est combine un encodage/decodage ambisonique ; 
- la figure 14 represente schematiquement le traitement 

15 des signaux issus d'un microphone comportant une pluralite 
de capteurs de pression agenc£s sur une sphere, a titre 
illustratif, par encodage ambisonique, egalisation et 
compensation de champ proche au sens de 1 1 invention. ' 

20 On se refere tout d'abord a la figure 1 qui represente a 
titre illustratif un systeme global de spatialisation 
sonore. Un module la de simulation d'une scene virtuelle 
definit un objet sonore comme une source virtuelle d'un 
signal, par exemple monophonique, de position choisie dans 

25 l'espace tridimensional et qui definit une direction du 
son. II peut etre prevu en outre des specifications de la 
geometrie d'une salle virtuelle, pour simuler une 
reverberation du son. Un module de traitement 11 applique 
une gestion d'une ou plusieurs de ces sources par rapport 

30 a un auditeur (definition d'une position virtuelle des 
sources par rapport a cet auditeur) . II met en oeuvre un 
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processeur d f effet de salle pour simuler des 
reverberations ou autres en appliquant des retards et/ou 
des filtrages usuels. Les signaux ainsi construits sont 
transmis a un module 2a d'encodage spatial des 
5 contributions elementaires des sources. 

Parallelement, une prise de son naturelle peut §tre 
effectuSe dans le cadre d f un enregistrement sonore par un 
ou plusieurs . microphones disposes de fagon choisie par 

10 rapport aux sources reelles (module lb) . Les signaux 
capt^s par les microphones sont encodes par un module 2b. 
Les signaux acquis et encodes peuvent §tre trans formes 
selon un format de representation intermediaire (module 
3b), avant d'etre mix<§s par le module 3 aux signaux 

15 generes par le module la et encodes par le module 2a (issu 
des sources virtuelles) . Les signaux mix#s sont ensuite 
transmis, ou encore memorises sur un support, en vue d'une 
restitution ulterieure (flSche TR) . lis sont ensuite 
appliques a un module de decodage 5, en vue de la 

20 restitution sur un dispositif de restitution 6 comportant 
des hauts-parleurs. Le cas 6ch6ant, l'Stape de decodage 5 
peut §tre preced^e d'une etape de manipulation du champ 
sonore, par exemple par rotation, grace a un module de 
traitement 4 prevu en amont du module de decodage 5 . 

25 

Le dispositif de restitution peut se presenter sous la 
forme d'une multiplicity de hauts-parleurs, agenc£s par 
exemple a la surface d'une sphere dans une configuration 
tridimensionnelle (p6riphonique) pour assurer, a la 
30 restitution, notamment un ressenti d'une direction du son 
dans I'espace tridimensional . A cet effet, un auditeur se 
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place generalement au centre de la sphere formee par le 
reseau de haut-parleurs, ce centre correspondant au point 
de perception auditive cite ci-avant. En variante, les 
hauts-parleurs du dispositif de restitution peuvent etre 
agences dans un plan (configuration panoramique 
bidimensionnelle) , les hauts-parleurs etant disposes en 
particulier sur un cercle et l'auditeur se placant 
habituellement au centre de ce cercle. Dans une autre 
variante, le dispositif de restitution peut se presenter 
sous la forme d'un dispositif de type "surround" (5.1). 
Enfin, dans une variante avantageuse, le dispositif de 
restitution peut se presenter sous la forme d'un casque a 
deux ecouteurs pour une synthese binaurale du son 
restitue, qui permet a l'auditeur de ressentir une 
15 direction des sources dans l'espace tridimensionnel, comme 
on le verra plus loin de facon detaillee. Un tel 
dispositif de restitution a deux hauts-parleurs, pour un 
ressenti dans l'espace tridimensionnel, peut se presenter 
aussi sous la forme d'un dispositif de restitution 
transaurale, a deux hauts-parleurs disposes a une distance 
choisie d'un auditeur. 



20 



On se ref£re maintenant a la figure 2 pour decrire un 
encodage spatial et un decodage pour une restitution 
25 sonore tridimensionnelle, de sources sonores elementaires . 
On transmet a un module d' encodage spatial 2 le signal 
issu d'une source 1 a N, ainsi que sa position (reelle ou 
virtuelle) . Sa position peut etre aussi bien definie en 
terme d' incidence (direction de la source vue de 
l'auditeur) qu'en terme de distance entre cette source et 
un auditeur. La plurality des signaux ainsi encodes permet 
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d'obtenir une representation multi-canale d'un champ 
sonore global. Les signaux encodes sont transmis (fleche 
TR) a un dispositif de restitution sonore 6, pour une 
restitution sonore dans 1 ■ espace tridimensional , comme 
5 indique ci-avant en reference a la figure 1. 

On se ref£re maintenant a la figure 3 pour deer ire ci- 
apres la representation ambisonique par des harmoniques 
sph£riques dans l 1 espace tridimensional, d'un champ 

10 acoustique. On considere une zone autour d f une origine O 
(sphere de rayon R) exempte de source acoustique. On 
adopte un systdme de coordonnSes spheriques dans lequel 
chaque vecteur r dds 1" origine O a un point de la sphere 
est d§crit par un azimut G r# une elevation 6 r et un rayon r 

15 (correspondant a la distance a 1' origine O) . 

Le champ de pression p(r) a 1 • int^rieur de cette sphdre 
(r < R oH R est le rayon de la sphere) peut s'ecrire dans 
le domaine frequent iel comme une serie dont les termes 

20 sont les produits ponder£s de f one t ions angulaires y% n (0,S) 
et de fonction radiale j m (kr) qui dependent ainsi d'un 
terme de propagation oil k=27tf/c / oil f est la frequence 
sonore et c est la vitesse du son dans le milieu de 
propagation. 

25 

Le champ de pression s 1 exprime alors par : 

Pir)=f J j m j m {kr) 2 Bg n Y^ D \e ri S r ) [Al] 
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L 1 ensemble des facteurs de ponderation B° m , qui sont 
implicitement fonction de la frequence, decrivent ainsi le 
champ de pression dans la zone consider#e. Pour cette 
raison, ces facteurs sont appeles "composantes harmoniques 
5 spheriques" et repr<§sentent une expression f requentielle 
du son (ou du champ de pression) dans la base des 

harmoniques spheriques Y% n . 

Les fonctions angulaires sont appelees "harmoniques 
10 spheriques" et sont d^finies par : 



[A2] 



V [m + n)\ [sin n0 si <7 = -1 

15 

ou 

PmnfsinS) sont des fonctions de Legendre de degre m et 
d'ordre n ; 

8 p , q est le symbole de Kronecker (€gal a 1 si p=q et 0, 
20 sinon) 

Les harmoniques spheriques forment une base orthonormee ou 
les produits scalaires entre composantes harmoniques et, 
de fa<?on genirale entre deux fonctions F et G, sont 
25 respectivement definies par : 

( Y mn \ Y £n' ^ = S mn? 5 nri 5 atf . [A 1 2] 
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Les harmoniques sph^riques sont des fonctions reelles 
bornees, comme represents sur la figure 4, en fonction de 
5 I 1 ordre m et des indices n et a. Les parties sombres et 
claires correspondent respectivement aux valeurs positives 
et negatives des fonctions harmoniques sph^riques . Plus 
1' ordre m est elevS et plus la frequence angulaire (et 
done la discrimination entre fonctions) est elevSe. Les 
10 fonctions radiales j m (kr) sont des fonctions de Bessel 
spheriques, dont le module est illustre pour quelques 
valeurs de 1' ordre m dans la figure 5. 

On peut donner une interpretation de la representation 
15 ambisonique par une base d 1 harmoniques spheriques comme 
suit. Les composantes ambisoniques de m§me ordre m 
expriment finalement des "dSrivees" ou des "moments" 
d 1 ordre m du champ de pression au voisinage de I'origine O 
(centre de la sphere representee sur la figure 3) . 

20 

En particulier, Bqo =W decrit la grandeur scalaire de la 

pression, tandis que B$ = X, Bff = Y, B$ = Z sont lies aux 
gradients de pression (ou encore a la velocity 
particulaire), a I'origine O. Ces quatre premieres 
15 composantes W, X, Y et Z sont obtenues lors d'une prise de 
son naturelle a l'aide de microphones omnidirectif s (pour 
la composante W d' ordre 0) et bidirectifs (pour les trois 
autres composantes suivantes) . En utilisant un plus grand 
nombre de transducteurs acoustiques, un traitement 
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d'avantage de composantes 
eleves superieurs a 1) . 
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egalisation, permet d'obtenir 
ambisoniques (ordres m plus 
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En prenant en compte des composantes supplementaires 
d' ordre plus eleve (superieur a 1) , done en augment ant la 
resolution angulaire de la description ambisonique, on 
accede a une approximation du champ de pression sur un 
voisinage plus large au regard de la longueur d'onde de 
l'onde sonore, autour de l'origine O. On comprendra ainsi 
qu'il existe une relation etroite entre la resolution 
angulaire (ordre des harmoniques spheriques) et la portee 
radiale (rayon r) qui peut etre representee. En bref, 
lorsque l'on s'ecarte spatialement du point d'origine O de 
la figure 3, plus le nombre de composantes ambisoniques 
est eleve (ordre M eleve) et meilleure est la 
representation du son par 1- ensemble de ces composantes 
ambisoniques. On comprendra aussi que la representation 
ambisonique du son est toutefois moins satisf aisante au 
fur et a mesure que l'on s'eloigne de l'origine O. Cet 
effet devient critique en particulier pour des frequences 
sonores elevees (de longueur d'onde courte) . On a done 
interet a obtenir un nombre de composantes ambisoniques 
qui soit le plus grand possible, ce qui permet de creer 
une region de l'espace autour du point de perception, dans 
laquelle la restitution du son est fidele et dont les 
dimensions sont croissantes avec le nombre total de 
composantes . 



30 



On decrit ci-apres une application a systeme 
d'encodage/ transmission/ restitution d'un son spatialise. 
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En pratique, un systeme ambisonique prend en compte un 
sous-ensemble de composantes harmoniques spheriques, comme 
d^crit ci-avant. On parle d'un systeme d'ordre M lorsque 
celui-ci prend en compte. des composantes ambisoniques 
d'indice m < M. Lorsqu'il s'agit d'une restitution par un 
dispositif de restitution a hauts-parleurs, on comprendra 
que si ces hauts-parleurs sont disposes dans un plan 
horizontal, seules les harmoniques d'indice m=n sont 
exploitees. En revanche, lorsque le dispositif de 
restitution comporte des hauts-parleurs disposes sur la 
surface d'une sphere ("pSriphonie") , on peut en principe 
exploiter autant d 1 harmoniques qu'il existe de haut- 
parleurs . 

On d^signe par la reference S le signal de pression porte 
par vine onde plane et captg au point O correspondant au 
centre de la sphere de la figure 3 (origine de la base en 
coordonnees spheriques). L 1 incidence de l'onde est decrite 
par l'azimut 9 et 1" elevation 5. L 1 expression des 
composantes du champ associ£ a cette onde plane est donnee 
par la relation : 



BZ„=SY£„(0,S) [A3] 

Pour encoder (simuler) une source en champ proche a une 

distance p de 1 'origine O, on applique un filtre i^ p/c) 
pour "incurver" la forme des fronts d'onde, en consid^rant 
qu'un champ proche 6met, en premiere approximation, une 
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onde spherique. Les composantes encodees du champ 
deviennent : 

5 

et 1« expression du filtre precit6 F^ pld) est donnee par la 
relation : 

^(?n-n)\n\ 

10 ou co = 27if est la pulsation de I'onde, f £tant la 
frequence du son. 

Ces deux dernieres relations [A4] et [A5] montrent 
finalement que, aussi bien pour une source virtuelle 

15 (simulee) que pour une source r€elle en champ proche, les 
composantes du son dans la representation ambisonique 
s'expriment math6matiquement (en particulier 

analytiquement) sous la forme d'un polynome, ici de 
Bessel, de puissance m et dont la variable (c/2jcop) est 

20 inversement proportionnelle a la frequence sonore. 



Ainsi, on comprendra que : 

- dans le cas d'une onde plane, I'encodage produit des 
signaux qui ne different du signal d'origine que d f un gain 

25 reel, fini, ce qui correspond a un encodage purement 
directionnel (relation [A3] ) ; 

- dans le cas d'une onde spherique (source en champ 
proche) , le filtre suppl^mentaire F^ p/c \co) encode 
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1 » information de distance en introduisant , dans 
1' expression des composantes ambisoniques, des rapports 
d« amplitudes complexes qui dependent de la frequence, 
comme exprime dans la relation [A5] . 

II est a noter que ce filtre supplementaire est de type 
"Integra teur", avec un effet d' amplification croissant et 
divergent (non-borne) au fur et a mesure que les 
frequences sonores decroissent vers zero. La figure 6 
montre, pour chaque ordre m, une augmentation du gain en 
basses frequences (ici la premiere distance p = i m ) . II 
s'agit done de filtres instables et divergents lorsque 
l'on cherche a les appliquer a des signaux audio 
quelconques. Cette divergence est d'autant plus critique 
15 pour les ordres m de valeur elevee. 

On comprendra en particulier, a partir des relations [A3] , 
[A4] , et [A5] , que la modelisation d'une source virtuelle 
en champ proche presente des composantes ambisoniques 

20 divergentes en basses frequences, de fagon 
particulierement critique pour des ordres m eleves, comme 
represents sur la figure 6. Cette divergence, dans les 
basses frequences, correspond au phenomene de "bass boost" 
enonce ci-avant. II se manifeste aussi en acquisition 

25 sonore, pour des sources reelles. 



Pour cette raison notamment, l'approche ambisonique, en 
particulier pour des ordres m eleves, n>a pas connu, dans 
l'etat de la technique, une application concrete (autre 
30 que theorique) dans le traitement du son. 
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On comprend en particulier qu'une compensation du champ 
proche est nScessaire pour respecter, a la restitution, la 
forme des fronts d'ondes encodes dans la representation 
ambisonique. En se ref grant a la figure 7, un dispositif 
5 de restitution comporte vine pluralite de hauts-parleurs 
HPi, disposes a une meme distance R, dans I'exemple 
decrit, d'un point de perception auditive P. Sur cette 
figure 7 : 

- chaque point oil se situe un haut-parleur HP ± correspond 
10 a un point de restitution enonce ci-avant, 

- le point P est le point de perception auditive 6nonc6 
ci-avant, 

- ces points sont separes de la seconde distance R enoncee 
ci-avant, 

15 tandis que sur la figure 3 decrite ci-avant : 

- le point 0 correspond au point de reference, Enonce ci- 
avant, qui forme I'origine de la base des harmoniques 
sph^riques , 

- le point M correspond a la position d'une source (r€elle 
20 ou virtuelle) situee a la premiere distance p, enoncee ci- 
avant, du point de reference O. 



Selon 1' invention, on introduit une pr£- compensation du 
champ proche au stade meme de l'encodage, cette 
25 compensation mettant en jeu des filtres de la forme 

analytique —J ^ et qui s'appliquent aux composantes 

F}h RIC \o>) 

ambisoniques B% n precitees. 
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Selon l ! un des avantages que procure 1 1 invention, 

1' amplification F^ ,c \a)) dont 1 1 ef f et apparalt sur la 
figure 6 est compensee par 1 1 attenuation du filtre 

applique des l'encodage 1 . En particulier, les 

F< R,c \a>) 

coefficients de ce filtre de compensation . } sont 

F!h Rlc \co) 

croissants avec la frequence du son et, en particulier, 
tendent vers z<§ro, pour les basses frequences. 
Avantageusement, cette pre- compensation, effectuee des 
l'encodage, assure que les donnees transmises ne sont pas 
divergentes pour les basses frequences. 

Pour indiquer la signification physique de la distance R 
qui intervient dans le filtre de compensation, on 
considere, a titre illustratif, une onde plane r£elle, 
initiale, k 1 1 acquisition des signaux sonores. Pour 
simuler un effet de champ proche de cette source 
lointaine, on applique le premier filtre de la relation 
[A5] , comme indique dans la relation [A4] . La distance p 
represente alors une distance entre une source virtuelle 
proche M et le point O representant l'origine de la base 
spherique de la figure 3. On applique ainsi un premier 
filtre de simulation de champ proche pour simuler la 
presence d'une source virtuelle a la distance p decrite 
ci-avant. NSanmoins, d'une part, comme indique ci-avant, 
les termes du coefficient de ce filtre divergent dans les 
basses frequences (figure 6) et, d 1 autre part, la distance 
p precit£e ne reprSsentera pas forcement la distance entre 
les hauts-parleurs d'un dispositif de restitution et un 
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point P de perception (figure 7). Selon 1 ■ invention, on 
applique une prS-compensation, a l'encodage, mettant en 

jeu un filtre de type 1 comme indiqug ci-avant, ce 

Ffn R,c \a>) 

qui permet, d'une part, de transmettre des signaux borngs, 
5 et, d' autre part, de choisir la distance R, d£s 
l'encodage, pour la restitution du son a partir des hauts- 
parleurs HPi, tel que represents sur la figure 7. En 
particulier, on comprendra que si l'on a simulS, a 
1" acquisition, une source virtuelle placee a la distance p 

10 de l'origine O, a la restitution (figure 7), un auditeur 
place au point P de perception auditive (§l une distance R 
des hauts-parleurs HP ± ) ressentira, a 1' audition, la 
presence d'une source sonore S, placSe a la distance p du 
point de perception P et qui correspond a la source 

15 virtuelle simulee lors de 1 1 acquisition. 



Ainsi, la pr£- compensation du champ proche des hauts- 
parleurs (places a la distance R) , au stade de l'encodage, 
peut etre combinee & un effet de champ proche simule d'une 
source virtuelle placee a une distance p. A l'encodage, on 
met finalement en jeu un filtre total resultant, d'une 
part, de la simulation du champ proche, et # d 1 autre part, 
de la compensation du champ proche, les coefficients de ce 
filtre pouvant s'exprimer analytiquement par la relation : 



Ffh R,c \co) 



[All] 
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Le filtre total donne par la relation [All] est stable et 
constitue la partie "encodage de distance" dans 1 ' encodage 
ambisonique spatial selon 1' invention, tel que represents 
sur la figure 8. Les coefficients de ces filtres 
correspondent a des fonctions de transfert monotones de la 
frequence, qui tendent vers la valeur 1 en • hautes 
frequences et vers la valeur (R/p) m en basses frequences. 
En se referant a la figure 9, les spectres d'energie des 

filtres H™ c (P ,c > r '°\g>) traduisent 1 • amplification des 
composantes encodees, dues a l'effet de champ de la source 
virtuelle (placee ici a une distance p = l m) , avec une 
pre -compensation du champ des hauts-parleurs (places a une 
distance R = 1,5 m) . L ' amplification en decibels est done 
positive lorsque p < R (cas de la figure 9) et negative 
quand p > R (cas de la figure 10 ou p = 3 m et R = l, 5 m) . 
Dans un dispositif de restitution spatialisee, la distance 
R entre un point de perception auditive et les haut- 
parleurs HPi est ef fectivement de l'ordre de un ou 
quel que s metres . 



En se referant a nouveau a la figure 8, on comprendra que, 
outre les parametres de direction G et 8 habituels, on 
transmettra une information sur les distances qui 
interviennent a 1- encodage. Ainsi, les fonctions 

25 angulaires correspondant aux harmoniques spheriques Y^ 6 ' 5 ^ 
sont conservees pour 1' encodage directionnel . 

Toutefois, au sens de la presente invention, on prevoit en 
outre des filtres totaux (compensation de champ proche et, 
30 le cas echeant, simulation d'un champ proche) 
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H NFC(p/c,R/c)^ gui sont appliques aux composantes 
ambisoniques, en fonction de leur ordre m, pour rSaliser 
l 1 encodage de la distance, comme represents sur la figure 
8. Un mode de realisation de ces filtres dans le domaine 
5 audionumerique sera decrit en detail plus loin. 

On remarquera. en particulier que ces filtres peuvent etre 
appliques des meme l 1 encodage de distance (r) et avant 
meme 1 1 encodage de direction (0,5). On comprendra ainsi 

10 que les etapes a) et b) ci -avant peuvent etre rassemblees 
en une meme etape globale, ou meme §tre interverties (avec 
un encodage de distance et filtrage de compensation, 
suivis d'un encodage de direction) . Le procede selon 
l 1 invention ne se limite done pas une mise en oeuvre 

15 successive dans le temps des etapes a) et b) . 

La figure 11A repr^sente une visualisation (vue de dessus) 
d'une reconstruction d'un champ proche avec compensation, 
d'une onde spherique, dans le plan horizontal (avec les 

20 m§mes paramdtres de distance que ceux de la figure 9) , 
pour un systdme d 1 ordre total M = 15 et une restitution 
sur 32 hauts-parleurs. Sur la figure 11B, on a represents 
la propagation de I'onde sonore initiale §l partir d'une 
source en champ proche situSe a une distance p d'un point 

25 de l'espace d' acquisition qui correspond, dans l'espace de 
restitution, au point P de la figure 7 de perception 
auditive. On remarque sur la figure 11A que les auditeurs 
(symbolisms par des tetes schematisees) peuvent localiser 
la source virtuelle en un meme lieu g£ographique situg a 

30 la distance p du point de perception P sur la figure 11B. 
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On verifie bien ainsi que la forme du front d'onde encode 
est respectee apres d£codage et restitution. Toutefois, on 
constate sensiblement des interferences k droite du point 
P tel que represents sur la figure 11A qui sont dues au 
5 fait que le nombre de hauts-parleurs (done de composantes 
ambisoniques prises en compte) n'est pas suffisant pour 
restituer parf aitement le front d f ondes en jeu sur toute 
la surface dSlimitee par les haut-parleurs . 

10 Dans ce qui suit, on d£crit, a titre d'exemple, 
l'obtention d r un filtre audionumerique pour la mise en 
ceuvre du procede au sens de 1 ! invention. 

Comme indique ci-avant, si l'on cherche a simuler un effet 
15 de champ proche, compens€ des l'encodage, on applique aux 
composantes ambisoniques du son un filtre de la forme : 



H NFC{plc > R/c) {(o) = 




[All] 



De 1' expression de la simulation d'un champ proche donne 
par la relation [A5] , il apparait que pour des sources 
lointaines (p = oo) , la relation [All] devient simplement : 



m '(») 

25 



F$«"ha» m ( } [A12] 



II apparait done de cette' derniere relation [A12] que le 
cas ou la source k simuler emet en champ lointain (source 
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lointaine) n'est qu'un cas particulier de 1 ' expression 
generale du filtre formulee dans la relation [All] . 

Dans le domaine des traitements audionumeriques, un 
procede avantageux pour definir un filtre numerique a 
partir de 1' expression analytique de ce filtre dans le 
domaine analogique a temps continu cons is te en une 
" trans formee bi 1 inSaire " . 

On exprime d'abord la relation [A5] sous la forme d'une 
transformed de Laplace, ce qui correspond a : 



Fm iP) -£j^W. (2tP) [A13] 

15 ou x = p/c (c etant la vitesse acoustique dans le milieu, 
typiquement 340 m/s dans l'air). 

La transformee bilineaire consiste a presenter, pour une 
frequence d- echantillonnage f s , la relation [All] sous la 
20 forme : 



[A14] 



si m est impair et 

25 
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si m est pair, 



1-z 



-1 



ou z est defini par p = 2f s - — par rapport a la relation 

1 + 2 1 

[A13] precedente, 



et avec 



et 



*0 =1-2 



Re (^ m 9 q) t \ X m,q 



a 



a 



x x = -2 



1- J 



m 9 g\ 



a 



« a 2 



10 



, 0 (^)/2 = l_fM et JC («+l)/2 



15 



oil a = 4f s R/c pour x=a 
et a = 4f s p/c pour x=b 

X m , q sont les q racines successives du polynome de Bessel 



m 



m 
q=l 
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et sont exprimees dans le tableau 1 ci-aprds, pour 
differents ordres m, sous les formes respectives de leur 
partie reelle, leur module (scares par une virgule) et 
leur valeur (reelle) lorsque m est impair. 
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Tableau 1 : valeurs R e [X m>q ] , \x m>g \ (et R e [X m ^ lorsque m 
est impair) d'un polynSme de Bessel calculees a l'aide du 
logiciel de calcul MATLAB©. 



m-1 


" -2.0000000000 ' 1 


m— 2 


" -3.0000000000, 3.4641016151 — 


m-3 


" -3.6778146454, 5.0830828022 ; -4.6443707093 " 


ixi — 4 


4.2075787944, 6.7787315854 ; -5.7924212056, 6.0465298776 


m-5 


4.6493486064, 8.5220456027 ; -6.7039127983, 7.5557873219 ; 
-7.2934771907 


m=6 


5.0318644956, 10.2983543043 ; -7.471416712 7, 9.1329783045 ; 
-8.4967187917, 8.6720541026 


m=7 


-S. 3713537579, 12.0990553610 
-9.5165810563, 10.1324122997 


? -8.1402783273, 10.7585400670 ; 
f -9.9435737171 


m-8 


-5.6779678978, 13.9186233016 
-10 .4096815813, 11.6507064310 


• —8.7365784344, 12.4208298072 ; 
; -11.1757720865, 11.3096817388 


m-9 


D.ybBb^l5964, 15.7532774523 ; -9.2768797744, 14.1121936859 ; 
-11.2088436390, 13.2131216226 ; -12.2587358086, 12.7419414392 ; 
-12.5940383634 


m-10 


6.2178324673, 17.6003068759 ; -9.77243 91337, 15.8272658299 ; 
-11.9350566572, 14.8106929213 ; -13.2305819310, 14.2242555605 ; 
-13 .8440898109, 13 . 9524261065 


m=ll | 


6.4594441798, 19.4576958063 ; -10.2312965678, 17.5621095176 ; 
-12.6026749098, 16.4371594915 ; -14.1157847751, 15.7463731900 ; 
-14.9684597220, 15.3663558234 ; -15.2446796908 


m-12 


6.6860466156, 21.3239012076 ; -10.6594171817, 19.3137363168 ; 
-13.2220085001, 18.0879209819 ; -14.9311424804, 17.3012295772 ; 
-15.9945411996, 16.8242165032 ; -16.5068440226, 16.5978151615 


m=13 


6.8997344413, 23.1977134580 ; -11.0613619668, 21.0798161546 ; 
-13.8007456514, 19.7594692366 ; -15.6887605582, 18.8836767359 ; 
-16.9411835315, 18.3181073534 ; -17.6605041890, 17.9988179873 ; 
-17.8954193236 


m-14 


7.1021737668, 25.0781652657 ; -11.4407047669, 22.8584924996 ; 
-14.3447919297, 21.4490520815 ; -16.3976939224, 20.4898067617 ; 
-17.8220011429, 19.8423306934 ; -18.7262916698, 19.4389130000 ; 
-19.1663428016, 19.2447495545 
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m-15 


-7.2947137247, 26.9644699653 ; -11.8003034312, 24.6482592959 ; ~ 

-14.8587939669, 23.1544615283 ; -17 . 0649i81370 , 22.1165594535 ; 
-18.6471986915, 21.3925954403 ; -19.7191341042, 20.9118275261 ; 
-20.3418287818, 20.6361378957 ; -20.5462183256 


m=16 


-7.4784635949, 28.8559784487 ;-12 . 1424827551 , 26.4478760957 ; 

-15.3464816324, 24.8738935490 ; -17.6959363478, 23.7614799683 ; 
-19.4246523327, 22.9655586516 ; -20.6502404436, 22.4128776078 ; 
-21.4379698156, 22.0627133056 ; -21.8237730778, 21.8926662470 


m-17 


7.6543475694, 30.7521483222 ; -12.4691619784, 28.2563077987 ; 
-15.8108990691, 26.6058519104 ; -18.2951775164, 25.4225585034 ; 
-20.1605894729, 24.5585534450 ; -21.5282660840, 23.9384287933 ; 
-22.4668764601, 23.5193877036 ; -23.0161527444, 23.2766166711 ; 
-23.1970582109 


m-18 


7.8231445835, 32.6525213363 ; -12.7819455282, 30.0726807554 ; 
-16.2545681590, 28.3490792784 ; -18.8662638563, 27.0981271991 ; 
-20.8600257104, 26.1693913642 ; -22.3600808236, 25.4856138632 ; 
-23.4378933084, 25.0022244227 ; -24.1362741870, 24.6925542646 ; 
-24.4798038436, 24.5412441597 


m=19 


7.9855178345, 34.5567065132 ; -13.0821901901, 31.8962504142 ; 

-16.6796008200, 30.1025072510 ; -19.4122071436, 28.7867778706 ; 
-21.5270719955, 27.7962699865 ; -23.1512112785, 27.0520753105 ; 
-24.3584393996, 26.5081174988 ; -25.1941793616, 26.1363057951 ; 
-25.6855663388, 25.9191817486 ; -25.8480312755 



On implemente ainsi les filtres numeriques, a partir des 
valeurs du tableau 1, en prevoyant des cascades de 
cellules d'ordre 2 (pour m pair), et une cellule 
supplementaire (pour m impair) , a partir des relations 
[A14] donnees ci-avant. 

On realise ainsi des filtres numeriques sous une forme de 
reponse impulsionnelle infinie, aisement parametrable 
comme montre ci-avant. II est a noter qu'une 
implementation sous une forme de reponse impulsionnelle 
finie peut etre envisagee et consiste a calculer le 
spectre complexe de la fonction de transfert a partir de 
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5 



15 



25 



la formule analytique, puis a en deduire une reponse 
impulsionnelle finie par transformee de Fourier inverse. 
On applique ensuite une operation de convolution pour le 
f iltrage . 



Ainsi, en introduisant cette pre-compensation du champ 
proche a I'encodage, on definit une representation 
ambisonique modifiee (figure 8) , en adoptant comme 
representation transmissible des signaux exprimes dans le 
10 domaine frequentiel, sous la forme : 



Comme indique ci-avant, R est une distance de reference a 
laquelle est associe un effet de champ proche compense et 
c est la vitesse du son (typiquement 340 m/s dans l'air) . 
Cette representation ambisonique modif i6e possdde les 
memes proprietes de scalabilite (schematiquement 
representee par des donnees transmises "entourees" pres de 
20 la fleche TR de la figure 1) et obeit aux memes 
transformations de rotation du champ (module 4 de la 
figure 1) que la representation ambisonique habituelle. 



On indique ci-apres les operations a mettre en ceuvre pour 
le decodage des signaux ambisoniques regus . 



On indique tout d'abord que 1" operation de decodage est 
adaptable a un dispositif de restitution quelconque, de 
rayon R 2/ different de la distance de reference R ci- 
30 avant. A cet effet, on applique des filtres de type 
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H NFC(p/c,R/c) 

(co) , tels que decrits plus haut, mais avec des 
parametres de distance R et R 2 , au lieu de p et R. En 
particulier, il est a noter que seul le parametre R/c est 
a m§moriser (et/ou transmettre) entre l'encodage et le 
5 decodage . 

En se rSferant a la figure 12, le module de filtrage qui y 
est represents est prSvu par exemple dans une unite de 
traitement d'un dispositif de restitution. Les composantes 
ambisoniques regues ont ete pre-compensees a l'encodage 
pour une distance de reference R x en tant que seconde 
distance. Toutefois, le dispositif de restitution comporte 
une pluralite de haut s -par leurs disposes a une troisieme 
distance R 2 d'un point de perception auditive P, cette 
troisieme distance R 2 etant differente de la seconde 
distance pr§cit§e R x . Le module de filtrage de la figure 

12 , sous la forme (co) , adapte alors , a la 

reception des donnSes, la pre- compensation a la distance 
Ri pour une restitution & la distance R 2 . Bien entendu, 
comme indique ci-avant, le dispositif de restitution 
regoit aussi le parametre Ri/c . 

Il est k noter que 1' invention permet en outre de mixer 
plusieurs representations ambisoniques de champs sonores 
25 (sources reelles et/ou virtuelles) , dont les distances de 
reference R sont diffSrentes (le cas echeant avec des 
distances de reference infinies et correspondant a des 
sources lointaines) . Pref erentiellement , on filtrera une 
pr€- compensation de toutes ces sources a une distance de 
30 reference la plus petite, avant de m^langer les signaux 
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ambisoniques, ce qui permet a la restitution d'obtenir une 
definition correcte du relief sonore. 

Dans le cadre d'un traitement dit de " focal i sat ion sonore" 
avec, a la restitution, un effet d' enrichissement sonore 
pour une direction choisie de l'espace (a la maniere d'un 
projecteur lumineux eclairant dans une direction choisie 
en optique) , impliquant un traitement matriciel de 
focalisation sonore (avec ponderation des composantes 
ambisoniques), on applique avantageusement l'encodage de 
distance avec pre -compensation de champ proche de facon 
combinee au traitement de focalisation. 

Dans ce qui suit, on decrit un procede de decodage 
ambi sonique, avec compensation du champ proche des hauts- 
parleurs, a la restitution. 

Pour reconstruire un champ acoustique encode suivant le 

formalisme ambisonique, a partir des composantes B%„ et 
en utilisant des hauts-parleurs d'un dispositif de 
restitution qui prevoit un emplacement "ideal" d'un 
auditeur qui correspond au point de restitution P de la 
figure 7, l'onde emise par chaque haut-parleur est definie 
par un traitement prealable de "rS-encodage" du champ 
ambisonique au centre du dispositif de restitution, comme 
suit. 

Dans ce contexte de "rS-encodage" , on considere dans un 
premier temps et pour simplification que les sources 
30 emettent en champ lointain. 
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En se referant a nouveau a la figure 7, l'onde emise par 
un haut-parleur d'indice i et d' incidence (9i et 5i) est 
alimente par un signal Si. Ce haut-parleur participe a la 

reconstruction de la composante B mn , par sa contribution 

Le vecteur ci des coefficients d'encodage associes aux 
hauts-parleurs d'indice i s'exprime par la relation : 



10 



Ci = 



[Bl] 



Le vecteur S des signaux emanant de 1' ensemble des N 
hauts-parleurs est donne par 1' expression : 



15 



S = 



'Si ' 

s 2 

_S N 



[B2] 



20 



La matrice d'encodage de ces N hauts-parleurs (qui 
correspond finalement a une matrice de "rS-encodage") , 
s'exprime par la relation : 



C=[C! C 2 .„C N ] 



[B3] 
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ou chaque terme Ci represente un vecteur selon la relation 
[Bl] ci-avant . 

Ainsi, la reconstruction du champ arabisonique B 1 est 
5 definie par la relation : 



ZJt+1 



0 mn 



= as 



[B4] 



La relation [B4] definit ainsi vine operation de re- 
10 encodage, pr€alable a la restitution. Finalement, le 
decodage, en tant que tel, consiste a comparer les signaux 
ambisoniques originaux et recus par le dispositif de 
• restitution, sous la forme : 



15 



-#11 



B = 



B 



11 



[B5] 



'mn 



aux. signaux re-encodes B, pour definir la relation 
general e : 

B' = B [B6] 
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II s'agit, en particulier, de determiner les coefficients 
d'une matrice de decodage D, qui verifie la relation : 



10 



S = D.B [B7] 

De preference, le nombre de hauts-parleurs est super ieur 
ou egal au nombre de composantes ambisoniques a decoder et 
la matrice de decodage D s'exprime, en fonction de la 
matrice de rg-encodage C, sous la forme : 

D = C T .(C.C T y l [B8 ] 



ou la notation C T correspond a la transposee de la matrice 
C. 

15 

II est a noter que la definition d ! un decodage v^rifiant 
des criteres diffgrents par bandes de frequences est 
possible, ce qui permet d'offrir une restitution optimis^e 
en fonction des conditions d'ecoute, notamment pour ce qui 
20 concerne la contrainte de positionnement au centre O de la 
sphere de la figure 3, lors de la restitution. A cet 
effet, on prevoit avantageusement un filtrage simple, en 
egalisation f requentielle par paliers, a chaque composante 
ambi sonique . 

25 

Toutefois, pour obtenir une reconstruction d'une onde 
originellement encodee, il faut corriger l'hypothdse de 
champ lointain pour les hauts-parleurs, c'est-a-dire 
exprimer l 1 effet de leur champ proche dans la matrice de 
30 r£-encodage C ci-avant et inverser ce nouveau syst£me pour 
definir le decodeur. A cet effet, en supposant une 
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concentricity des hauts-parleurs (disposes a une meme 
distance R du point P de la figure 7) , tous les hauts- 
parleurs ont un meme effet de champ proche F m R/c ^(a>) , sur 
chaque composante ambisonique du type B'Z. 



mn 



En 



5 introduisant les termes de champ proche sous la forme 



d'une matrice diagonale, la relation 
devient : 



[B4] ci-avant 



10 



B>=Diag{\ if' <» Ff'ico) -F Rlc {co) F* ,c {a» ••]). 



La relation [B7] ci-avant devient : 



C.S [B9] 



S = DJDiag 



I Ff'ico) Ff'ia) F£' c (o>) F* ,c {<o) 



,R/c 



R/c, 



J J 



JB 



[BIO] 



15 



20 



25 



Ainsi, 1» operation de matricage est precedee par une 
operation de filtrage qui compense le champ proche sur 

chaque composante B%„ , et qui peut etre mise en oeuvre 
sous forme numerique, comme decrit ci-avant, en reference 
a la relation [A14] . 

On retiendra qu'en pratique, la matrice C de "re-encodage" 
est propre au dispositif de restitution. Ses coefficients 
peuvent itre determines initialement par parametrage et 
caracterisation sonore du dispositif de restitution 
reagissant a un excitation predeterminee . La matrice de 
decodage D est, elle aussi, propre au dispositif de 
restitution. Ses coefficients peuvent etre determines par 
la relation [B8] . En reprenant la notation precedente ou 
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B est la matrice des composantes ambisoniques pre- 
compensees, ces dernieres peuvent etre transmises au 
dispositif de restitution sous forme matricielle B avec : 



B = Diag 



.B 



Le dispositif de restitution decode ensuite les donnees 
recues sous forme matricielle B (vecteur colonne des 
composantes transmises) en appliquant la matrice de 
decodage D aux composantes ambisoniques pre-compensees, 
pour former les signaux Si destines a alimenter les haut- 
parleurs HPi, avec : 



15 



20 



S = 



Si 



= D.B 



[Bll] 



En se referant a nouveau a la figure 12, si une operation 
de decodage doit etre adaptee a un dispositif de 
restitution de rayon R 2 different de la distance de 
reference r 1# un module d- adaptation prealable au decodage 
proprement dit et decrit ci-avant permet de filtrer chaque 



composante ambisonique B%„ , pour 1' adapter a un 
dispositif de restitution de rayon R 2 . L' operation de 
decodage proprement dite est effectuee ensuite, comme 
decrit ci-avant, en reference a la relation [Bll] . 
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On decrit ci-apres une application de 1 » invention a la 
synthase binaurale. 

On se ref£re k la figure 13A sur laquelle un auditeur 
5 disposant d'un casque a deux ecouteurs d'un dispositif de 
synthase binaurale est represents. Les deux oreilles de 
l 1 auditeur sont disposees a des points respectifs O l 
(oreille gauche) et Or (oreille droite) de l'espace. Le 
centre de la tete de 1' auditeur est dispose au point O et 

10 le rayon de la t£te de 1' auditeur est de valeur a. Une 
source sonore doit etre pergue audit ivement a un point M 
de l'espace, situe a une distance r du centre de la t§te 
de 1' auditeur (et respectivement a des distances r R de 
l 1 oreille droite et r L de 1' oreille gauche). Par ailleurs, 

15 la direction de la source placee au point M est definie 
par les vecteurs r , r R et r L . 

De faq:on generale, la synthese binaurale se definit comme 
suit. 

20 

Chaque auditeur a une forme d 1 oreille qui lui est propre. 
La perception d'un son dans 1 ' espace par cet auditeur se 
fait par apprentissage, depuis la naissance, en fonction 
de la forme des oreilles (notamment la forme des pavilions 

25 et les dimensions de la tete) propre a cet auditeur. La 
perception d'un son dans 1' espace se manifeste entre 
autres par le fait que le son parvient a une oreille, 
avant 1' autre oreille, ce qui se traduit par un retard x 
entre les signaux k emettre par chaque ecouteur du 

30 dispositif de restitution appliquant la synthese 
binaurale. 
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Le dispositif de restitution est parametre initialement , 
pour un m§me auditeur, en balayant une source sonore 
autour de sa tete, a une m§me distance R du centre de sa 
5 t§te. On comprendra ainsi que cette distance R peut etre 
consideree comme une distance entre un "point de 
restitution" comme enonce ci-avant et un point de 
perception auditive (ici le centre O de la t§te de 
1 1 auditeur) . 

10 

Dans ce qui suit, . 1 » indice L est associ€ au signal a 
restituer par l'£couteur accole k 1' oreille gauche et 
1' indice R est associe au signal a restituer par 
l'6couteur accol£ a l 1 oreille droite. En se ref grant a la 
15 figure 13B, on applique au signal initial S un retard pour 
chaque voie destinee & produire un signal pour un ecouteur 
distinct. Ces retards x L et x R sont fonction d f un retard 
maximum Xmax qui correspond ici au rapport a/c oil a, comme 
indique precedemment , correspond au rayon de la tete de 

20 l f auditeur et c a la vitesse du son. En particulier, ces 
retards sont dgfinis en fonction de la difference de 
distance du point O (centre de la tete) au point M 
(position de la source dont le son est a restituer, sur la 
figure 13A) et de chaque oreille a ce point M. 

25 Avantageusement, on applique en outre des gains respect if s 
gi, et g R/ a chaque voie, qui sont fonction d f un rapport 
des distances du point O au point M et de chaque oreille 
au point M. Des modules respectifs appliques a chaque voie 
2 L et 2 R encodent les signaux de chaque voie, dans une 

30 representation ambisonique, avec pr£- compensation de champ 
proche NFC (pour "Near Field Compensation") au sens de la 



WO 2004/049299 



CT/FR2003/003367 



48 



presente invention. On comprendra ainsi que, par la mise 
en ceuvre du procede au sens de la presente invention, on 
peut definir les signaux issus de la source M, non 
seulement par leur direction (angles azimutaux 9 L et 9 R et 
5 angles d- elevation 5 L et 8 R ) , mais aussi en fonction de la 
distance separant chaque oreille r L et r R de la source M. 
Les signaux ainsi encodes sont transmis au dispositif de 
restitution comportant des modules de decodage 
ambisonique, pour chaque voie, 5 L et 5 R . Ainsi, on 
10 applique un encodage/decodage ambisonique, avec 
compensation de champ proche, pour chaque voie (ecouteur 
gauche, ecouteur droit) dans la restitution avec synthese 
binaurale (ici de type "B-FORMAT") , sous forme dedoublee. 
La compensation de champ proche s'effectue, pour chaque 
voie, avec comme premiere distance p une distance r L et r R 
entre chaque oreille et la position M de la source sonore 
a restituer. 



15 



20 



25 



30 



On decrit ci-apres une application de la compensation au 
sens de 1' invention, au contexte de 1 ' acquisition sonore 
en representation ambisonique. 

On se refere a la figure 14 sur laquelle un microphone 141 
comporte une plurality de capsules transductrices , 
capables de capter des pressions acoustiques et restituer 
des signaux electriques S lt ...,s N . Les capsules CAPi sont 
agencees sur une sphere de rayon r predetermine (ici, une 
sphere rigide, telle qu'une balle de ping-pong par 
exemple) . Les capsules sont espacees d'un pas regulier sur 
la sphere. En pratique, on choisit le nombre N de capsules 
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en fonction de l'ordre M desire pour la . representation 
ambisonique . 

On indique ci-apres, dans le contexte d'un microphone 
5 comportant des capsules agencees sur une sphere rigide, 
comment compenser l'effet de champ proche, des 1 1 encodage 
dans le contexte ambisonique. On montrera ainsi que la 
pre -compensation du champ proche peut s'appliquer non 
seulement pour la simulation de source virtuelle, comme 
10 indiquS ci-avant, mais aussi a 1 1 acquisition et, de fagon 
plus generale, en cpmbinant la pre -compensation de champ 
proche a tous types de traitements impliquant une 
representation ambisonique . 

15 En presence d'une sphere rigide (susceptible d'introduire 
une diffraction des ondes sonores regues) , la relation 
[Al] donnee. ci-avant devient : 



^r(g/)-S - i^, : Z*™ 7 ™**,) [CI] 

tr=±l 

20 

Les d^rivees des fonctions de Hankel spheriques h" m 
obeissent a la loi de recurrence : 



(2m + l)h-(x) — m h m _i(x) — (m + 1) h m+ \{pc) [C2] 

25 



On deduit les composantes ambisoniques B% n du champ 
initial £ partir du champ de pression a la surface de la 
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sphdre, en mettant en oeuvre des operations de projection 
et d 1 Sgalisation donnees par la relation : 



10 



*Sn=EQm<PrT£n> 4 * [C3] 

Dans cette expression, EQ m est un filtre egaliseur qui 
compense une ponderation W m qui est li€e a la directivity 
des capsules et qui inclut en outre la diffraction par la 
sphere rigide. 

L 1 expression de ce filtre EQ m est donnee par la relation 
suivante : 



EQ m =^- = (kr) 2 h^(kr)j- m+1 [C4] 

15 

Les coefficients de ce filtre d ' egalisation ne sont pas 
stables et on obtient un gain infini en tres basses 
frequences- D'ailleurs, il convient de noter que les 
composantes harmoniques spheriques, elles-memes , ne sont 
20 pas d 1 amplitude finie lorsque le champ sonore n'est pas 
limite a une propagation d f ondes planes, c'est-a-dire 
issues de sources lointaines, comme on l'a vu 
precedemment . 



25 Par ailleurs, si, plutot que de prevoir des capsules 
encastrees dans une sphere solide, on prevoit des capsules 
de type cardioldes, avec une directivite en champ lointain 
donnee par l f expression : 
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G(0) = a + (l-a)cos0 [C5] 

En considerant ces capsules mont^es sur un support 
"transparent acoustiquement ", le terme de ponderation a 
5 compenser devient : 

W m =j m (a jm{kr)-j{\-a)jm\krj) [ C 6] 

II apparait encore que les coefficients d'un filtre 
10 d' egalisation correspondant §l 1 1 inverse analytique de 
cette ponderation donnee par la relation [C6] sont 
divergents pour les tres basses frequences, 

De fagon generale, on indique que pour tout type de 
15 directivite de capteurs, le gain du filtre EQ m pour 
compenser la ponderation W m liee a la directivite des 
capteurs est infini pour les basses frequences sonores. En 
se refer ant a la figure 14, on applique avantageusement 
une pr6- compensation de champ proche dans l 1 expression 
20 meme du filtre d' egalisation EQ m , donnee par la relation : 



25 



Fi R ' C \o>) 

Ainsi, les signaux S x k S N sont recupSr^s du microphone 
141. Le cas echeant, on applique une pre- egalisation de 
ces signaux par un module de traitement 142. Le module 143 
permet d'exprimer ces signaux dans le contexte 
ambisonique, sous forme matricielle. Le module 144 
applique le filtre de la relation [C7] aux composantes 
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ambisoniques exprim^es en fonction du rayon r de la sphdre 
du microphone 141. La compensation de champ proche 
s'effectue pour une distance de reference R en tant que 
seconde distance. Les signaux encodes et ainsi filtr^s par 
5 le module 144 peuvent §tre transmis, le cas 6ch6ant, avec 
le paramStre repr§sentatif de la distance de reference 
R/c. 

Ainsi, il apparalt dans les diff brents modes de 

10 realisation lies respectivement k la creation d'une source 
virtuelle en champ proche, a 1 ' acquisition de signaux 
sonores issues de sources rSelles, ou merae a la 
restitution (pour compenser un effet de champ proche des 
hauts-parleurs) , que la compensation de champ proche au 

15 sens de la pr^sente invention peut s'appliquer a tous 
types de traitements faisant intervenir vine representation 
ambisonique. Cette compensation de champ proche permet 
d'appliquer la representation ambisonique a une 
multiplicity de contextes sonores ou la direction d'une 

20 source et avantageusement sa distance doivent §tre prises 
en compte. De plus, la possibility de la representation de 
phenomenes sonores de tous types (champs proches ou 
lointains) dans le contexte ambisonique est assuree par 
cette pre -compensation, du fait de la limitation a des 

25 valeurs reelles finies des composantes ambisoniques. 

Bien entendu, la presente invention ne se limite pas a la 
forme de realisation d£crite ci-avant k titre d'exemple ; 
elle s'etend a d'autres variantes . 
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Ainsi, on comprendra que la prS-compensation de champ 
proche peut etre integree, a l'encodage, autant pour vine 
source proche que pour une source lointaine. Dans ce 
dernier cas (source lointaine et reception d'ondes 
5 planes) , la distance p exprimee ci-avant sera consideree 
comme infinie, sans modifier de fa?on substantielle 
1 'expression des filtres H ra donnee ci-avant. Ainsi, le 
traitement utilisant des processeurs d'effet de salle qui 
fournissent en general des signaux decorrelSs utilisables 
10 pour modeliser le champ diffus tardif (reverberation 
tardive) peut etre combine a une pre -compensation de champ 
proche. On peut considerer que ces signaux sont de meme 
Snergie et correspondent a une part de champ diffus 

correspondant a la composante omnidirective W = Bqq 
15 (figure 4) . On peut alors construire les di verses 
composantes harmoniques spheriques (avec un ordre M 
choisi) en appliquant une correction de gain pour chaque 
composante ambisonique et on applique une compensation de 
champ proche des hauts-parleurs (avec une distance de 
20 reference R separant les haut-parleurs du point de 
perception auditive comme represents sur la figure 7) . 

Bien entendu, le principe d'encodage au sens de la 
prSsente invention est generalisable a des modules de 
25 rayonnement autres que des sources monopolaires (reelles 
ou virtuelles) et/ou des hauts-parleurs. En effet, toute 
forme de rayonnement (notamment une source Stalee dans 
I'espace) peut §tre exprim<§e par integration d'une 
distribution continue de sources elementaires ponctuelles. 
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En outre, dans le contexte de la restitution, il est 
possible d' adapter la compensation de champ proche a tout 
contexte de restitution. A cet effet, il peut etre prevu 
de calculer des fonctions de transfert (re-encodage des 
composantes harmoniques spheriques de champ proche pour 
chaque haut-parleur, compte tenu d'une propagation reelle 
dans la salle ou le son est restitue) , ainsi qu'une 
inversion de ce re-encodage pour redefinir le decodage. 

On a deer it ci-avant un procede de decodage dans lequel on 
appliquait un systeme matriciel faisant intervenir les 
composantes ambisoniques . Dans une variante, il peut etre 
prevu un traitement generalise par transformees de Fourier 
rapides (circulaire ou spherique) pour limiter les temps 
de calcul et les ressources informatiques (en terme de 
memoire) necessaires au traitement de decodage. 

Comme indique ci-avant en reference aux figures 9 et 10, 
on constate que le choix d'une distance de reference R par 
rapport a la distance p de la source en champ proche 
introduit une difference de gain pour differentes valeurs 
de la frequence sonore. On indique que le procede 
d'encodage avec pre-compensation peut etre couple a une 
compression audionumerique permettant de quantifier et 
d'ajuster le gain pour chaque sous-bande frequentielle. 

Avantageusement, la presente invention s' applique a tous 
types de systemes de spatialisation sonore, notamment pour 
des applications de type "rSalitS virtuelle" (navigation 
dans des scenes virtuelles dans l'espace tridimensionnel, 
jeux avec spatialisation sonore tridimensionnelle, 
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conversations de type "chat" sonorisees sur le reseau 
Internet), a des sonif ications d» interfaces, a des 
logiciels d' edition audio pour enregistrer, mixer et 
restituer de la musique, mais aussi a 1 1 acquisition, a 
5 partir d 1 usage de microphones tridimensionals, pour la 
prise de son musicale ou cinematographique, ou encore pour 
la transmission d 1 ambiance sonore sur Internet, par 
exemple pour des "WebCam" sonorisees. 
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RE VEND I CAT IONS 

1. Procede de traitement de donnees sonores, dans lequel : 

a) on code des signaux represent at if s d'au moins un son se 
5 propageant dans l'espace tridimensionnel et issu d'une 

source situee a une premiere distance (p) d'un point de 
reference (O) , pour obtenir une representation du son par 
des composantes (B. 9 ) exprimees dans une base 
d'harmoniques spheriques, d'origine correspondant audit 
10 point de reference (o) , 

b) et on applique auxdites composantes (B^) une 
compensation d'un effet de champ proche par un filtrage 
qui est fonction d'une seconde distance (R) definissant 
sensiblement, pour une restitution du son par un 
dispositif de restitution, une distance entre un point de 
restitution (HPJ et un point (P) de perception auditive. 

2. Procede selon la revendication 1, dans lequel, ladite 
source etant lointaine du point de reference (o) , 

- on obtient des composantes d'ordres successifs m pour la 
representation du son dans ladite base d'harmoniques 
spheriques, et 

- on applique un f litre (l/FJ dont les coefficients, 
appliques chacun a une composante d'ordre m, s'expriment 
analytiquement sous la forme de 1 ' inverse d'un polyndme de 
puissance m, dont la variable est inversement 
proportionnelle a la frequence sonore et a ladite seconde 
distance (R) , pour compenser un effet de champ proche au 
niveau du dispositif de restitution. 
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3. Proc€de selon la revendication 1, dans lequel # ladite 
source £tant vine source virtuelle pr§vue 3. ladite premiere 
distance (p) , 

- on obtient des composantes d'ordres successifs m pour la 
5 representation du son dans ladite base d 1 harmoniques 

spheriques , et 

- on applique un filtre global (H m ) dont les coefficients, 
appliques chacun a une composante d'ordre m, s'expriment 
analytiquement sous la forme d'une fraction, dont : 

10 - le numerateur est un polyndme de puissance m, dont 

la variable est inversement proportionnelle a la 
frequence sonore et a ladite premiere distance (p) , 
pour simuler un effet de champ proche de la source 
virtuelle, et 

15 ~ 1© denominateur est un polynome de puissance m, 

dont la variable est inversement proportionnelle a la 
frequence sonore et a ladite seconde distance (R) , 
pour compenser 1' effet du champ proche de la source 
virtuelle dans les basses frequences sonores. 

20 

4. Precede selon l'une des revendications precedentes, 
dans lequel on transmet au dispositif de restitution les 
donn6es cod^es et filtr£es aux etapes a) et b) avec un 
parametre representatif de ladite seconde distance (R/c) . 

25 

5. Proced£ selon l'une des revendications 1 a 3, dans 
lequel, le dispositif de restitution comportant des moyens 
de lecture d'un support memoir e, on memorise sur un 
support m€moire destine a gtre lu par le dispositif de 

30 restitution les donn£es codees et filtr^es aux Stapes a) 
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et b) avec un parametre reprSsentatif de ladite seconde 
distance (R/c) • 

6. Procede selon l'une des revendi cat ions 4 et 5, dans 
lequel, prealablement a une restitution sonore par un 
dispositif de restitution comportant une plurality de 
hauts-parleurs disposes a une troisieme distance (R 2 ) 
dudit point de perception auditive (P) , on applique aux 
donnges cod^es et filtrees un filtre d' adaptation 
(H m (R1/c ' R2/c) ) dont les coefficients sont fonction desdites 
seconde (R x ) et troisieme distances (R 2 ) . 

7. Procede selon la revendication 6, dans lequel les 
coefficients dudit filtre d' adaptation ( Hm <Ri/c,R2/c) j ^ 
appliques chacun a une composante d'ordre m, s'expriment 
analytiquement sous la forme d'une fraction, dont : 

" le numerateur est un polynome de puissance m, dont la 
variable est inversement proportionnelle a la frequence 
sonore et a ladite seconde distance (R) , 

~ et le denominateur est un polynome de puissance m, dont 
la variable est inversement proportionnelle a la frequence 
sonore et a ladite troisieme distance (R 2 ) . 

8. Procede selon l'une des revendi cat ions 2, 3 et 7, dans 
lequel, pour la mise en ceuvre de l'etape b) , on prevoit : 

- pour des composantes d'ordre m pair, des filtres 
audionumeriques sous la forme d'une cascade de cellules 
d'ordre deux ; et 

- pour des composantes d'ordre m impair, des filtres 
audionumeriques sous la forme d 1 une cascade de cellules 
d'ordre deux et une cellule suppl€mentaire d'ordre un. 
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9. Proced£ selon la revendication 8, dans lequel les 
coefficients d»un filtre audionumerique, pour une 
composante d'ordre m, sont definis a partir des valeurs 

5 num^riques des racines desdits polyndmes de puissance m. 

10. Proced€ selon l'une des revendications 2, 3, 7, 8 et 
9 # dans lequel lesdits polynomes sont des polynomes de 
Bessel . 
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11. Proc£de selon l'une des revendications 1, 2 et 4 a 10, 
dans lequel on prevoit un microphone comportant un reseau 
de transducteurs acoustiques agenc^s sensiblement sur la 
surface d'une sphere dont le centre correspond 
sensiblement audit point de reference (O) , pour obtenir 
lesdits signaux representatif s d'au moins un son se 
propageant dans l'espace tridimensional. 

12. Proc^de selon la revendication 11, dans lequel on 
applique a l'etape b) un filtre global pour, d'une part, 
compenser un effet de champ proche en fonction de ladite 
seconde distance (R) et, d' autre part, €galiser les 
signaux issus des transducteurs pour compenser une 
ponderation de directivity desdits transducteurs. 

13. Procede selon l'une des revendications 11 et 12, dans 
lequel on prevoit un nombre de transducteurs fonction d'un 
nombre total choisi de composantes pour repr£senter le son 
dans ladite base d' harmoniques spheriques. 
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14. Proced£ selon 1'une des revendications precedentes, 
dans lequel on choisit a l'gtape a) un nombre total de 
composantes dans la base des harmoniques sph€riques pour 
obtenir, a la restitution, une region de 1 1 espace autour 
du point de perception (P) dans laquelle la restitution du 
son est fidele et dont les dimensions sont croissantes 
avec le nombre total de composantes. 

15. Procede selon la revendication 14, dans lequel on 
prevoit un dispositif de restitution comportant un nombre 
de haut-parleurs au moins egal audit nombre total de 
composantes . 



16. Procede selon l'une des revendications 1 a 5 et 8 a 
15 13, dans lequel : 

- on prevoit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a une 
distance choisie d'un auditeur, 

- on obtient, pour cet auditeur, une information de 
20 ressenti de la position dans l'espace de sources sonores 

situees a une distance de reference predetermine (R) de 
1" auditeur, et 

- on applique la compensation de l'etape b) avec ladite 
distance de reference sensiblement en tant que seconde 

25 distance. 

17. Proced# selon l'une des revendications 1 & 3 et 8 a 
13, prises en combinaison avec l'une des revendications 4 
et 5, dans lequel : 
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- on pr^voit un dispositif de restitution comportant au 
moins un premier et un second haut-parleur disposes a une 
distance choisie d'un auditeur 7 

- on obtient, pour cet auditeur, une information de 
5 ressenti de la position dans l'espace de sources sonores 

situ£es a une distance de reference predeterminee (R 2 ) de 
1 1 auditeur, et 

- prealablement a une restitution sonore par le dispositif 
de restitution, on applique aux donn^es cod^es et filtr€es 

10 aux etapes a) et b) un filtre d' adaptation (H m (R/c ' R2/c) ) 
dont les coefficients sont fonction de la seconde distance 
(R) et sensiblement de la distance de reference (R 2 ) . 

18. Proced£ selon l'une des revendications 16 et 17, dans 
15 lequel : 

- le dispositif de restitution comporte un casque & deux 
ecouteurs pour les oreilles respectives de 1' auditeur, et 

- separement pour chaque gcouteur, on applique le codage 
et le filtrage des etapes a) et b) pour des signaux 

20 respectifs destines a alimenter chaque ecouteur, avec, en 
tant que premiere distance (p) , respectivement une 
distance (r R ,r L ) separant chaque oreille d'une position 
(M) d'une source a restituer. 

25 19. Procede selon l'une des revendications precedentes, 
dans lequel on met en forme, aux etapes a) et b) , un 
systeme matriciel comportant au moins : 

- une matrice (B) comportant lesdites composantes dans la 
base des harmoniques sph6riques, et 
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- une matrice diagonale (Diag(l/FJ) dont les coefficients 
correspondent a des coefficients de filtrage de l ! 6tape 
b), 

et on multiplie lesdites matrices pour obtenir une matrice 
5 r£sultat de composantes compensees (B). 



2 0. Proc^de selon la revendication 19, dans lequel : 

- le dispositif de restitution comporte une pluralite de 
haut-parleurs disposes sensiblement & une meme distance 

10 (R) du point de perception auditive (P) , et 

- pour decoder lesdites donnees codees et filtrees aux 
etapes a) et b) et former des signaux adapt^s pour 
alimenter lesdits haut-parleurs : 

* on forme un systeme matriciel comportant ladite 
15 matrice resultat (B) et une matrice de d<§codage (D) 

predeterminee, propre au dispositif de restitution, 
et 

* on obtient une matrice (S) comportant des 
coefficients representatif s des signaux 

20 d 1 alimentation des hauts-parleurs par multiplication 

de la matrice des composantes compensees (B) par 
ladite matrice de d£codage (D) . 



21. Dispositif d 1 acquisition sonore, comportant un 
microphone muni d'un reseau de transducteurs acoustiques 
disposes sensiblement sur la surface d'une sphdre, 
caract£ris€ en ce qu'il comporte en outre une unit£ de 
traitement agencee pour : 

- recevoir des signaux emanant chacun d'un transducteur, 
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- appliquer auxdits signaux un codage pour obtenir vine 
representation du son par des composantes (3^°) exprimees 
dans une base d'harmoniques sph€riques, d'origine 
correspondant au centre de ladite sphere (O) , 

- et appliquer auxdites composantes (B^ 0 ) un filtrage qui 
est fonctipn, d'une part, d'une distance correspondant au 
rayon de la sphere (r) et, d' autre part, d'une distance de 
reference (R) . 

22. Dispositif selon la revendication 21, caracteris6 en 
ce que ledit filtrage consiste, d'une part, a egaliser, en 
fonction du rayon de la sphere, les signaux issus des 
transducteurs pour compenser une pond^ration de 
directivite desdits transducteurs et, d' autre part, a 
compenser un effet de champ proche en fonction d'une 
distance de reference choisie (R) , def inissant 
sensiblement , pour une restitution du son, une distance 
entre un point de restitution (HP ± ) et un point (P) de 
perception auditive. 
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